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We utilize classical molecular dynamics to study the the quality (Q)-factors of monolayer CVD- 
grown graphene nanoresonators. In particular, we focus on the effects of intrinsic grain boundaries 
of different orientations, which result from the CVD growth process, on the Q-factors. For a range of 
misorientations orientation angles that are consistent with those seen experimentally in CVD-grown 
graphene, i.e. 0° to '~-^20°, we find that the Q-factors for graphene with intrinsic grain boundaries 
are 1-2 orders of magnitude smaller than that of pristine monolayer graphene. We find that the 
Q-factor degradation is strongly influenced by both the symmetry and structure of the 5-7 defect 
pairs that occur at the grain boundary. Because of this, we also demonstrate that find the Q-factors 
CVD-grown graphene can be significantly elevated, and approach that of pristine graphene, through 
application of modest (1%) tensile strain. 



I. INTRODUCTION 

Since its recent discovery as the simplest two- 
dimensional crystal structure^, graphene has been exten- 
sively studied not only for its unusual physical proper- 
ties resulting from its two-dimensional structure^Ti^, but 
also for its potential as the basic building block of fu- 
ture applications, i.e. nanoclcctromechanical systems 
(NEMS)^-ia. 

Graphene is viewed as an ideal material for NEMS- 
based sensing and detection applications due to its com- 
bination of extremely low mass and exceptional mechan- 
ical propertiesii; we note the recent review of Barton et 
al^ in this regard. However, one key issue limiting the 
applicability of graphene as a sensing component is its 
low quality (Q)-factor; the Q-factors of a 20-nm thick 
multilayer graphene sheet were found to range from 100 
to 1800 as the temperature decreased from 300 K to 
50 K^. Similarly low Q-factors between 2 and 30 were 
also observed by Sanchez et al^ for multilayer graphene 
sheets, while higher Q-factors with values up to 4000 
were reported using multilayer graphene oxide films^. 
Theoretically, the Q-factors of graphene were recently 
studied using classical molecular dynamics (MD) simu- 
lationai^^, where spurious edge modes that arc present 
in suspended graphene were proposed to have a key role 
in the low Q-factors that were observed experimentally. 
This hypothesis was recently validated by Barton et ali^, 
who found Q-factors approaching 2000 for graphene res- 
onators that were clamped on all slides, thus eliminating 
the spurious edge modes. 

These early experimental works on graphene nanores- 
onators utilized graphene flakes made via the scotch- 
tape method, which produces graphene sheets of vary- 
ing thickness, though importantly, each graphene layer 
in the sheet is single crystalline. However, graphene re- 
search has been transformed by the recent development 
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of the chemical-vapor-deposition (CVD) growth process 
to synthesize large area graphene sheetsi^"— . CVD tech- 
nology has the potential to revolutionize graphene-based 
sensing technology due to the resulting promise of wafer- 
scale graphene devices comprised of large arrays of single- 
layer graphene resonatorsi^. However, while large area 
graphene films are desirable for these and other graphene- 
based applications^^, the graphene films that result from 
CVD growth are polycrystalline, and thus are composed 
of many interconnected single crystalline graphene grains 
that intersect at grain boundaries having a range of mis- 
orientation angles2i, which generally are less than 20°—. 
While these grain boundaries are sometimes viewed fa- 
vorably for tunable electronic devices^^, they are likely 
to have a deleterious effect for graphene nanoresonators 
because the misorientation and the resulting non-ideal 
bonding at the grain boundary causes an increase in 
phonon scattering, which creates another energy dissi- 
pation mechanism and a lower Q-factor. 

Some very recent experimental studies, such as those of 
van der Zandc et alJ^ and Barton et alt^ have studied 
CVD-grown graphene nanoresonators, and have found 
relatively high Q-factors on the order of about 2000. 
However, those works also improved the Q-factors by re- 
moving spurious edge modesi^, and therefore the intrin- 
sic losses that occur in CVD-grown graphene due to the 
existence of the grain boundaries is unknown. Further- 
more, previous theoretical studiesi^ii^i^ focused on the 
energy dissipation mechanisms in pure graphene without 
grain boundaries. Therefore, the objective of the present 
work is to quantify, via classical MD simulations, the in- 
trinsic dissipation mechanisms introduced in CVD-grown 
graphene nanoresonators due to the presence of the grain 
boundaries. 



II. SIMULATION METHODOLOGY 

The starting point of our simulations is to note that re- 
cent experimental studies have found that the misorienta- 
tion angles at grain boundaries in CVD-grown graphene 
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FIG. 1: (Color online) Schematic of the graphene 
monolayers containing grain boundaries with various 
misorientation angles of (a) 0°;(b)1.4°; (c) 5.31°; (d) 

9.83°; (e) 12.83°; (f) 16.62°. 5-7 defect pairs are 
highlighted in yellow. 



lie mostly between 0° to 20°—. To systematically study 
the effects of different grain boundary orientations on the 
intrinsic loss mechanisms in graphene monolayers, we cre- 
ated graphene monolayers with a single grain boundary 
along the center line of the monolayer that runs along the 
armchair orientation with 6 different misorientation an- 
gles of 0° (pristine graphene), 1.4°, 5.31°, 9.83°, 12.83° 
and 16.62°. These configurations are shown in Fig. [1] 
where the diameter of all the graphene monolayers was 4 
nm. The initial configuration including the single grain 
boundary was generated by rotating two semicircular 
graphene monolayers and then piecing them together. At 
that point, energy minimization using the conjugate gra- 
dient algorithm was employed to optimize the structure, 
where the carbon-carbon interactions were modeled us- 
ing the AIREBO potential of Stuart et al.— , which is 
able to accurately simulate the forming and breaking of 
carbon bonds. A key point to note is that the creation of 
the grain boundary results in the formation of 5-7 defect 
ring pairs, which are colored in yellow in Fig. [TJ We can 
observe that as the grain boundary misorientation angle 
increases, the density of the 5-7 ring pairs increases. For 
example, for the smallest misorientation angle of 1.4° in 
Fig. [ija), there exists only a single 5-7 unit ring pair 
near the center of the monolayer. In contrast, as the 
grain boundary misorientation angle increases, the den- 
sity of the 5-7 ring pairs increases, as observed for the 
other cases in Fig. [T] 

The Q-factors were calculated and the intrinsic energy 
dissipation mechanisms studied using classical MD via 
the publicly available simulation code LAMMPS^. After 
obtaining the equilibrium graphene monolayer structures 
with the various grain boundary misorientation angles as 
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FIG. 2: Out-of-planc buckling of 9.83° defective 
graphene sheet after thermal equilibration. 
Visualization performed using VMD^. 



shown in Fig. [U we then performed a thermal equilibra- 
tion using a Nose-Hoover thermostat^l for 500 ps using 
a time step of 1 fs, i.e. within an NVT ensemble. During 
the equilibration, the edges of the graphene sheet were 
constrained in plane while the rest of the sheet was left 
free to move. Previous theoretical^^ and experimental^^ 
works have demonstrated that spurious edge vibrational 
modes, which arise due to the undcrcoordinatcd nature 
of bonding at the edges of graphene, have a dominant role 
in reducing the Q-factors of suspended graphene nanorcs- 
onators; these edge modes would also be present for sus- 
pended CVD-grown nanoresonators, and thus should be 
eliminated to maximize the Q-factor. Because of this, af- 
ter the thermal equilibration, the edges of the graphene 
monolayer were clamped at the equilibrium diameter that 
is established during the thermal equilibration to elimi- 
nate the possibility of spurious edge modes. 

During the thermal equilibration, out-of-plane buck- 
ling was observed due to the 5-7 defect pairs as observed 
in Fig. [2] such buckles were previously observed in MD 
simulations by Liu et al.— . We find that for the small- 
est grain boundary angle of 1.4° at 3K, the height of the 
buckle is about 3.3A. For larger grain boundary misori- 
entation angles, there are more buckles along the grain 
boundary due to the larger number of 5-7 defect pairs, 
which interact and lead to a decrease in the buckling 
height. The impact of the buckles on the Q-factors will 
be elucidated later; in particular, we will demonstrate 
the utility of tensile mechanical strain in enhancing the 
Q-factor by flattening out the buckles. 

After the thermal equilibration, the graphene mono- 
layer was actuated by assigning an initial sinusoidal ve- 
locity profile that ranged from zero at the clamped edges 
to a maximum at the center of the circular monolayer, 
and where the initial velocity was applied only in the 
vertical z-direction to be perpendicular to the graphene 
sheets, and was chosen to be sufficiently small such that 
the resulting oscillation of the graphene monolayer would 
be purely harmonic, i.e. the resulting increase in total 
energy due to the applied sinusoidal velocity was only 
about 0.1%. While the buckling results in a non-planar 
graphene monolayer, for consistency, the direction of the 
applied initial velocities were the same for both pristine 
and defective graphene sheets. After the velocity pro- 
file was prescribed, the resulting free oscillation of the 
graphene monolayer was performed within an energy con- 
serving (NVE) ensemble for 3000 ps. 

The Q-factors were calculated following the procedure 
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described by Vallabhaneni et al.— and Chu et al.— . 
Specifically, as described by Vallabhaneni et al.— , the 
variation of the displacement of the center of mass of the 
graphene monolayer was tracked for the duration of the 
simulation after the initial velocity is applied. The decay 
in the root mean square displacement was then fit to the 
following exponential curve {e~'^"*), which is then related 
to Q via Q= 0.5/(^, where ^ — 7/0; is the damping ratio 
and to is the angular vibrational frequency. Specifically, 
because the vibrational motion is predominately in the z- 
direction in our MD simulations, with little contribution 
from the motion in the x and y directions, we used the 
center of mass in the ^;-direction only to fit the damping 
curve. This approach is utilized in the present work as it 
avoids the necessity of extracting the external energy as 
previously performed by Kim and Park— J^*^. Figs. [3] 
and |4] show typical center of mass and natural frequency 
results for pristine and 5.31° defective graphene sheets, 
respectively, that were used to obtain the Q-factors. 

We make three other relevant comments here. First, 
we chose to study graphene monolayers with a sin- 
gle grain boundary rather than study polycrystalline 
graphene with a distribution of grain boundary misorien- 
tation angles. By comparing the Q-factors of monolayer 
graphene with a single grain boundary and by varying 
the misorientation angle of the single grain boundary to 
pristine graphene, and by utilizing temperatures ranging 
from ~0K to 300K, we aim to quantify the effects of each 
grain boundary orientation on the intrinsic loss mecha- 
nisms, or Q-factor. Second, the models in this work rep- 
resent extreme cases in that the defects run through the 
entire diameter of the graphene monolayer. Finally, we 
will also study the effects of tensile mechanical strain in 
enhancing the Q-factors of graphene. Previous studies 
on both graphenei^i2i and other nanostructures^ have 
demonstrated the effectiveness of strain in increasing the 
Q-factors. We will demonstrate the utility of strain in 
mitigating the effects of out of plane buckling due to the 
5-7 unit ring defects along the grain boundaries, thus el- 
evating the Q-factors. 



III. NUMERICAL RESULTS AND DISCUSSION 

Before presenting the misorientation angle-dependent 
results for the Q-factors, we first note again that the 5- 
7 defects along the grain boundary induce out of plane 
buckling in the graphene monolayer as shown in Fig. [21 
More specifically, we first show that both the orientation 
of the buckle (i.e. up or down), and its location along 
the grain boundary, significantly impacts not only the 
structural symmetry of the graphene sheet, but also the 
Q-factors; this is true for all grain boundary misorienta- 
tion angles. 

To demonstrate this, we considered several possible 
spatial locations of the single 5-7 defect pair for the 1.4° 
misorientation angle, as illustrated in Fig. [5] The first 
two cases, i.e. those depicted in Figs. [5ja) and (b). 
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FIG. 3: Center of mass and natural frequency for 
pristine graphene at 3K. 
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FIG. 4: Center of mass and natural frequency for 
defective graphene with a grain boundary 
misorientation angle of 5.31° at 3K. 



result in a significant asymmetry of the actuated oscilla- 
tion. However, the buckle in the third case in Fig. [SJc) 
is near the edge of the graphene monolayer, and thus 
better preserves the circular oscillation symmetry, and 
results in a Q-factor of about 10,000, which is almost 
1 order of magnitude larger than found for the cases in 
Fig. [5ja) and (b). However, to make the results of all 
the grain boundary misorientation cases to be compara- 
ble, we have placed the single 5-7 defect pair for the 1.4° 
case in the center of the monolayer as shown in Fig. [TJa) . 

To illustrate the infiuence of the orientation of the 
buckling, for each misorientation angle (except for the 
1.4° sheet with one buckle which thus has only two pos- 
sible buckling orientations, up or down), we tested five 
cases with different buckling patterns at 3K, where the 
different buckling patterns were generated by using differ- 
ent random velocity seeds during the thermal equilibra- 
tion. From the results shown in Fig. [51 we can see that 
all cases have a Q-factor that is 1-2 orders of magnitude 



4 



(al 



FIG. 5: (Color online) Schematic figures showing 
out-of-planc bucklings of 1.4° sheet with different defect 
pair positions after relaxation. 
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FIG. 6: (Color online) Variation of the Q-factor for 
monolayer graphene at 3K for different grain boundary 
misorientation angles. 



smaller than pristine, monolayer graphene (0°). Further- 
more, the range between the highest and lowest Q-factor 
for each misorientation angle spans approximately one 
order of magnitude. For consistency, we utilize the con- 
figuration with the highest Q-factor for each misorienta- 
tion angle, as indicated by the dashed line in Fig. [51 for 
the remainder of this study. 

The Q-factors of unstrained graphene, both pristine 
and with the different grain boundary orientations, are 
shown in Fig. [7] as a function of temperature. We find 
that the Q-factors of pristine graphene follow the rela- 
tionship Q ^ r~", where we find a = 1.28 for the present 
case, which is similar to the result of Kim and Park- 
when the different approaches for calculating Q are ac- 
counted for. We can see that for all temperatures, pris- 
tine graphene has a higher Q-factor than CVD graphene, 
where the difference in Q-factor at low temperatures is 
between 1-2 orders of magnitude, as shown in Fig. [6] We 
also observe that the Q factors for all graphene sheets 
with grain boundaries obey a similar functional form 
{Q ~ with respect to temperature, while for both 
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FIG. 7: (Color online) Variation of the Q-factor for 
unstrained monolayer graphene as a function of 
temperature and grain boundary misorientation angle. 



pristine and defected graphene sheets the Q factor drops 
rapidly by the time room temperature is reached. 

More interestingly, we find that the relationship be- 
tween the grain boundary misorientation angle and the 
Q-factor is non-monotonic. This is illustrated in Fig. |6l 
where the Q-factor is seen to first increase with increasing 
misorientation angle and then decrease when the angle 
becomes larger than about 10°. We believe this is due to 
a competition between two effects. 

First, as the grain boundary misorientation angle in- 
creases, the density of 5-7 defect pairs along the grain 
boundary increases, which should result in a decrease in 
the Q-factor with increasing grain boundary misorien- 
tation angle. On the other hand, while there are more 
buckles in the graphene monolayer with increasing grain 
boundary misorientation angle, we find that due to the 
interaction of the increasing number of defect pairs with 
increasing grain boundary misorientation angle that the 
circular symmetry of the graphene monolayer is better 
preserved. Specifically, as seen in Fig. [T] for the 1.4° 
case, a single, non-symmetric 5-7 defect pair exists near 
the center of the graphene sheet, which we have already 
discussed results in a substantial reduction in Q-factor. 
However, increasing the misorientation angle results in 
more 5-7 defects on each side of the central 5-7 defect, 
which enhances the overall symmetry of the graphene 
sheet. 

Furthermore, the height of the buckles tends to de- 
crease with increasing grain boundary misorientation an- 
gle, i.e. the buckling heights for the 1.4°, 5.31°, 9.83°, 
12.83° and 16.62° graphene sheets at 3K after relaxation 
are found to be 3.3A, 2.5A, 1.9A, 1.5Aand 2.2A, respec- 
tively. Thus, the increase in Q-factor due to the greater 
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FIG. 8: (Color online) Out-of-plane buckling of 16.62° 
defective graphene sheet after thermal equilibration. 



structural symmetry with increasing misorientation an- 
gle coupled with the corresponding reduction in buck- 
ling height counteracts the decrease in Q due to the in- 
crease in the defect density. To illustrate this concept, 
we show in Fig. [S]tlic equilibrium buckled configuration 
for the 16.62° case, where it can be seen that each of the 
four buckles forms from the combination of two smaller 
buckles. As a result, the buckles not only have larger 
buckling heights, but also break symmetry due to the 
fact that each buckle is composed of a smaller and larger 
sub-buckle. These factors couple to result in a smaller Q- 
factor for the 16.62° case, and the overall non-monotonic 
trend seen in Fig. [T] 

Due to the deleterious effect of the grain boundaries 
on the Q-factor, we also examine how the Q-factors of 
CVD-grown graphene can be enhanced. As suggested 
in various worksiSiHii^i^^ tensile mechanical strain is an 
effective approach to enhancing the Q-factors of nanos- 
tructures and NEMS. In our MD simulations, we imposed 
a modest, experimentally-accessible 1% tensile strain^i^ 
that was applied symmetrically (radially) outward from 
the center of the CVD-grown graphene sheets prior to 
the thermal equilibration and subsequent velocity-driven 
actuation. 

As shown by comparison of Figs. [7]and[9l the Q-factors 
of strained graphene with grain boundaries can, in some 
cases, approach those of pristine graphene. In addition, 
as shown in Fig. [TUl tensile strain increases the Q-factors 
of graphene with grain boundaries much more than pris- 
tine graphene, where the tensile strain-induced Q-factor 
enhancement for the graphene with grain boundaries can 
be larger than one order of magnitude. Fig. [S]also shows 
that strained graphene follows the same Q ^ T^" rela- 
tionship as unstrained graphene, though with a slightly 
lower exponent a of 1.11. 

We have found through our MD simulations that ten- 
sile strain increases the Q-factor by both increasing the 
natural frequency w, while simultaneously suppressing 
the damping 7 in the expression C. = j/ui, where Q= 
0.5/^. We found in analyzing the Q-factors that the ten- 
sile strain increased the natural frequency uj about the 
same amount for both pristine graphene and graphene 
with grain boundaries. However, the damping factor 
7 showed a significantly greater reduction for graphene 
with grain boundaries as compared to pristine graphene. 
This is because the tensile strain reduces the out of plane 
buckling that results due to the 5-7 defect pairs, thus 
further enhancing the structural integrity and suppress- 
ing the damping of the oscillating graphene monolayer. 
The specific buckling heights at 3K due to 1% tensile 
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FIG. 9: (Color online) Variation in the Q-factor for 
monolayer graphene under 1% tensile strain as a 
function of temperature and grain boundary 
misorientation angle. 
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FIG. 10: (Color online) Variation in the Q-factor as a 
function of grain boundary misorientation angle at 3K 
due to 1% tensile strain, and the the ratio the Q-factors 
under tensile strain to the Q-factors of unstrained 
graphene. 



strain for the grain boundary misorientation angles of 
1.4°, 5.31°, 9.83°, 12.83° and 16.62°: 2.5A, I.SA, 1.4A, 
l.lAand 1.7 A, respectively; these are clearly smaller than 
when no tensile strain is applied. 

We also note that when tensile strain is applied the 
misorientation angle- dependent Q-factor follows a simi- 
lar trend in Fig. [10] as was previously observed for the 
unstrained graphene in Fig. [71 For strained graphene at 
3K, we find that graphene with a misorientation angle 
of 12.83° has the highest Q factor, while the graphene 
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sheets with a 9.83° misorientation angle has a similarly 
high Q-factor. Overall, this demonstrates the important 
fact that CVD-grown graphene can, under modest, ex- 
perimentally accessible tensile strains, exhibit compara- 
ble performance for NEMS sensing applications as pris- 
tine graphene. 

IV. CONCLUSIONS 

We have utilized classical MD simulations to quantify 
the effects of grain boundaries in CVD-grown graphene 
on the Q-factors of graphene nanoresonators. Graphene 
with grain boundaries exhibit Q-factors that are 1-2 or- 
ders of magnitude smaller than pristine graphene. How- 
ever, the Q-factors follow a non-monotonic dependence 



on the grain boundary misorientation angle due to the 
competing effects of increased 5-7 defect pair density 
on one hand, and the increased structural symmetry 
and reduction in out of plane buckling heights on the 
other hand. Furthermore, for practical applications, the 
Q-factors of CVD-grown graphene can be enhanced by 
about one order of magnitude through the application 
of 1% tensile strain, which results in Q-factors that ap- 
proach those of pristine graphene. 
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